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Summary 


Spectral analysis and sensitivity analysis are used to discover which meteorological input parameters and road 
thermal properties are most important to a road surface temperature prediction model. The results show that the 
model is most sensitive to air temperature and cloud cover, and not very sensitive to road thermal properties. 


1. Introduction 

Numerical road ice prediction models have success- 
fully reproduced the major features of road surface 
conditions (Thornes 1984, 1989, Rayer 1987, Shao 
1990). However, few studies, if any, have looked at the 
interrelationships between road surface temperature 
and the geographical, road construction and meteor- 
ological inputs to the models, for instance as shown in 
Table I. How does the output of a numerical model 
respond to the model parameters and model input, 
whose values usually cover a considerable range? 
Spectral analysis and sensitivity analysis are used in this 
paper in order to answer some of these questions. 


Some general research has been carried out relating 
road surface temperatures to topography and certain 
weather parameters, such as the influence of topography 
and wind on the variation of minimum air temperatures 
and road surface temperatures (Tabony 1985, Bogren 
and Gustavsson 1988, Gustavsson 1990). Weather 
conditions have been classified for convenience into 
damped, intermediate and extreme days by Thornes 
(1989). Farmer and Tonkinson (1989) have looked at the 
sensitivity of the Meteorological Office model (Rayer 
1987) to road thermal properties and road surface 
albedo. These researches have shown the value of 





Table I. 


Parameters that control road surface temperature 





Geographical 


Road construction 


Meteorological 





Latitude 

Altitude 
Topography 

Sky view 

Land use 
Pollution 
Roughness length 
De-icing chemical 


Emissivity 
Albedo 
Traffic 


Depth of construction 
Thermal conductivity 
Thermal diffusivity 


Solar radiation 
Terrestrial radiation 
Air temperature 
Cloud cover 

Wind velocity 
Humidity 
Precipitation 
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investigating the interaction and sensitivity of road 
surface temperatures to their environment, but they lack 
sufficient systematic and quantitative information for 
model users, and especially for modellers who wish to 
improve their model predictions. This paper aims to 
examine the importance of the meteorological inputs 
and road thermal properties in order to further our 
understanding of the problems involved. 

A method of cross-spectral analysis is used to explore 
the cross-correlations between road surface temperatures 
and several meteorological parameters. Then a variety 
of sensitivity tests are carried out using the Icebreak 
model (Shao 1990) to examine the response of the model 
output to changes in both meteorological and road 
thermal properties, illustrating each with results from 
two control runs with typical values. Such results reflect 
model variability and provide suggestions to give better 
operational results, e.g. via a cautious selection of 
meteorological, road construction and geographical 
input data. 


2. Cross-correlation analysis 

As mentioned by Thornes and Shao (1991), model 
prediction error originates largely from meteorological 
input parameters, which include air temperature, dew- 
point, wind speed, cloud amount and type, and rainfall. 
All of these parameters are considered as random 
parameters or part of a time series. To examine 
correlations between road surface temperature and the 
meteorological parameters over different time-scales, 
spectral analysis is used. Details of the method are given 
in the Appendix. 

The results of the spectral analysis are listed in 
Tables II and III. There are 456 sets of hourly 
observations from Chapman’s Hill on the M5 motor- 
way including road and air temperatures, relative 
humidity and wind speed. The maximum lag considered 
in Table II is 24 hours. The number of cloud observ- 
ations is 129, and the maximum lag considered in 
Table III is also 24 hours. In these tables, T is the 
periodicity in hours; C’(7’) is the coherence squared for 
each time-period (coherence measures the degree of 


correlation between two time-series) and L(7) is the 
time difference of two time-series on the cycle of Tth 
periodicity in hours. The asterisk (*) in these tables 
indicates that the coherence passes the significance test 
at either a significance level of 0.01 (Table II) or 0.05 
(Table III). 

The results show that road surface temperature (T7;) 
has a close relationship with all the meteorological 
parameters except cloud cover on a half-day (12-hour) 
cycle. The small values in L(7T) indicate that the 
variation of 7, responds to that of meteorological 
parameters very quickly. A significant correlation also 
exists on some other cycles such as daily and 6-hour 
cycles. The results show that the air temperature is the 
most important factor to which road surface temperature 
responds. When there is a negative phase-difference the 
air temperature changes are driving the road surface 
temperatures (e.g. 6- and 16-hour lag), whereas with a 
positive value of phase difference the road surface 
temperature is driving air temperature changes (e.g. 8- 
and 24-hour lag). The lags are not simple to interpret in 
physical terms due to the effect of traffic, but basically 
one would expect that the lags should be positive during 
the day when the road is likely to be warmer than the air, 
and the opposite at night. 


3. Sensitivity tests for meteorological inputs 

The most reliable information about the physical 
nature of the variation of road surface temperature is 
obtained by actual measurement. Such an experimental 
investigation is expensive on time and resources and, in 
most cases, impossible. For instance one would have to 
wait a long time to measure all cloud combinations at 
a given wind speed. Sensitivity analysis, by using a 
mathematical and physical model, is a valuable way to 
quickly gain an insight into the controlling aspects of the 
variation of road surface temperature and to develop 
confidence in model predictions. The Icebreak model is 
run with standard meteorological observations (for 
the so-called control runs) for 4 December 1987 and 
24 February 1988, representing a ‘damped’ day and an 





Table Il. 


Chapman’s Hill, winter 1987/88, 456 samples 


Results of cross-power spectrum analysis between road 
surface temperature (X; (#)) and meteorological parameters (X; (t)) for 


Table Ill. Results of cross-power spectrum 
analysis between road surface temperature (X; (1) 
and cloud cover (X; (t)) for Chapman’s Hill, winter 





T 
hour 


Air temperature 


c’(T) 


L(T) 


Humidity 


c’(T) 


L(T) 


Wind speed 


c?(T) 


1987/88, 129 samples 





L(T) f Cloud amount Cloud type 





24.0 
16.0 
12.0 
8.0 
6.0 
4.0 
3.0 
2.0 


0.669* 
0.619* 
0.744* 
0.569* 
0.588* 
0.115 

0.147 

0.375* 


0.44 
—0.79 
0.08 
0.72 
—0.70 
—0.24 
—0.30 
0.02 


0.163 
0.096 
0.386* 
0.199 
0.267* 
0.020 
0.296* 
0.161 


—2.29 
2.19 
0.35 
0.99 

—0.76 

—0.22 

—0.24 
0.24 


0.067 
0.005 
0.277* 
0.021 
0.183 
0.129 
0.043 
0.129 


hour rahe | gl L(T) 


ca) L(T) 





1.73 
2.47 24.0 
0.38 20.0 
—1.16 15.0 
—0.50 12.0 
—0.23 10.0 
—0.33 8.0 
—0.05 6.0 


0.249* —0.22 
0.106 0.23 
0.210 0.83 
0.156 —0.34 
0.310* 0.57 
0.095 0.26 
0.554* 0.22 


0.197 0.49 
0.048 0.35 
0.013 —1.04 
0.057 —0.02 
0.048 —0.18 
0.132 0.62 
0.071 —0.06 





* see text for further explanation 


* see text for further explanation 
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‘extreme’ day. Figs. 1 and 2 show the results of model 
prediction with standard input. 

In the tests, each of the input parameters (air 
temperature, dew-point, wind speed, cloud amount and 
cloud type) is varied through a range of values and all 
other parameters are held constant. The range of values 
shown here for each parameter are added or subtracted 
for each 3-hour input. The model is run from noon to 
noon and input values are supplied for 1500, 1800, 2100, 
0000, 0300, 0600, 0900 and 1200rhrs: 

(a) Air temperature: + | °C: the results are shown in 

Figs 3 and 4. 

(b) Dew-point: + | °C (Figs 5 and 6). 

(c) Wind speed is the most variable atmospheric 

parameter considered. In the test, it undergoes a 

+ 3 kn change which is in accordance with typical 

forecast input errors (see Thornes and Shao (1991)). 

Figs 7 and 8 show the results. 

(d) Cloud amount: + | okta when the amount is less 

than 8 oktas (Figs 9 and 10). 

(e) Cloud type: +1 cloud type (1 =low cloud to 

3 = high cloud) where appropriate (Figs 11 and 12). 


The results of the test runs are displayed in Tables IV 
and V, and compared with the control runs. The bias, 
standard deviation (SD) and difference of minimum 
surface temperatures (D7 (min)) are given in the tables. 
The running difference shown in the tables is given by: 


running difference = Vn-V- 


where JV, is the bias, or SD, or DT (min) with modified 
input, and V. is that of control run. 

It is seen from the tables that of all of these 
meteorological parameters, air temperature is the most 


important. A change of + | °C in air temperature results 
in an extra error of about + 0.8 °C in the road surface 
temperature bias and minimum road surface temper- 
ature. This result accords with that of the spectral 
analysis discussed above. Cloud cover is the second 
most important factor. Its impact on the model 
predictions, especially on the prediction of minimum 
surface temperature, is as expected. The impact of wind 
speed is small and it is greater on an ‘extrem:=" day than 
on a ‘damped’ day. Dew-point seems to only have a 
minor influence on model output. 


4. Sensitivity tests for road thermal 
properties 

It is very difficult to obtain accurate values for the 
thermal properties of a road construction. Firstly, if 
sample cores are taken from a road to be analysed in the 
laboratory, inevitably water is needed to cool the 
diamond cutter. Water contamination in the sample is 
impossible to correct for, and significantly affects the 
thermal properties of the sample. Secondly, the depth of 
the sub-layers below the surface is often unknown, and 
therefore one has to make assumptions about the overall 
thermal properties of a road structure including the 
sub-base. Hence one has to take published data from the 
literature and select sensible values for the particular 
road being dealt with. Model runs are then used to test 
the sensitivity of road surface temperature to the 
thermal properties of the road construction. 

Table VI shows the typical values of road thermal 
properties used by several authors. The thermal 
properties of asphalt, concrete and soil are also unlikely 
to be the same at different sites, but inevitably the same 
values have to be used. It can be seen that in general the 





Table IV. Results of sensitivity tests for meteorological parameters and 
comparison to standard results for Chapman’s Hill, 4 December 1987 





Tests Output of tests 


Bias SD 


Running difference 


DT(min) A B © 





Standard 
T, (+1) 
T, (—1) 
DP (+1) 
DP (-—1) 
WS (+3) 0.07 0.278 
WS (—3) 0.24 0.313 
CA (+1) — — 
CA (-1) 0.04 0.466 
CT (+1) —0-14 0.432 
CT (-1) _ — 


0.15 
0.86 
—0.64 
0.15 
0.05 


0.289 
0.287 
0.266 
0.288 
0.246 


—0.81 


—0.16 


—0.27 
—0.48 


0.00 _ 
0.74 0.74 
—0.81 
0.03 
—0.16 
0.01 
0.04 


0.71 
—0.79 
0.00 
—0.10 
—0.08 
0.09 


0.03 


0.01 
0.04 
—0.11 —0.27 
-6.2 —0.48 





Key: 7, air temperature (°C) 
DP dew-point (°C) 
ws wind speed (kn) 
CA cloud amount in oktas 


cr cloud type (0 = no cloud, | = low cloud, 2 = medium cloud and 3 = high 


cloud 
A difference of mean errors 


B difference of standard deviations 
Cc difference of minimum surface temperature errors 
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Figure 1. Output with standard input for Chapman’s Hill 
4 December 1987. 
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Figure 3. Output with changes in air temperature for Chapman’s 
Hill 4 December 1987. 
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Figure 5. Output with changes in dew-point for Chapman’s Hill 
4 December 1987 


conductivity of asphalt is lower than that of concrete, 
and that heat capacity varies less than conductivity. 
For this sensitivity analysis the Icebreak model is run 
with a range of + 20% of the thermal parameters. The 
road construction is assumed to be: the top sub-layer 
(0-10cm asphalt), the middle sub-layer (10-40 cm 
concrete) and the bottom sub-layer (40-100 cm soil or 
granite). Table VII shows the values of conductivity and 
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Figure 2. Output with standard input for Chapman’s Hill 
24 February 1988. 
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Figure 4. Output with changes in air temperature for Chapman’s 
Hill 24 February 1988. 
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Figure 6. Output with changes in dew-point for Chapman’s Hill 
24 February 1988. 


capacity of asphalt, concrete and soil (or granite) used in 
the sensitivity tests. 

In Table VII, tests (a)-(b) and tests (c)-(d) are for a 
change of conductivity and capacity of asphalt and 
concrete, test (e) is for sandy soil and test (f) is for 
granite. The results of the tests are compared with those 
of the standard input in Tables VIII and IX. Comparing 
the model predictions of the control runs to those of the 
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Figure 7. Output with changes in wind speed for Chapman’s Hill 
4 December 1987. 
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Figure 9. Output with changes in cloud amount for Chapman’s Hill 
4 December 1987. 
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Figure 11. Output with changes in cloud type for Chapman’s Hill 
4 December 1987. 


modified runs, it can be seen that a 20% change of 
conductivity or capacity of the road sub-layers does not 
significantly influence model cutput. 


5. Conclusions 

Of the many factors that influence the variation of 
road surface temperature, air temperature is the most 
important. In routine operation of an ice prediction 
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Figure 8. Output with changes in wind speed for Chapman’s Hill 
24 February 1988. 
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Figure 10. Output with changes in cloud amount for Chapman’s Hill 
24 February 1988. 
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Figure 12. Output with changes in cloud type for Chapman's Hill 
24 February 1988. 


model when forecast meteorological input is provided, 
any error in the forecast of air temperature will result in 
a similar error in model road surface temperature 
prediction. 

Cloud cover (i.e. cloud amount and cloud type) is 
another key factor governing the variation in road 
surface temperatures. On an ‘extreme’ day, the road 
surface temperature responds quickly to a change of 
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Table V. Results of sensitivity tests for meteorological parameters and 
comparison to standard results for Chapman’s Hill, 24 February 1988 





Tests Output of tests 


SD 


Running difference 


Bias DT(min) A B © 





Standard 
T, (+1) 

st) 

DP (+1) 
DP (-—1) 
WS (+3) 
WS (3) 
CA (+1) 
CAC) 
CT (+1) 
CT (-—1)* 


—0.16 

0.68 
—0.99 
—0.16 
=i? 
—0.28 

0.02 
—0.28 
—0.09 
—0.18 
—0.06 


0.573 
0.571 
0.568 
0.571 
0.567 
0.638 
0.615 
0.738 
0.687 
0.620 
0.600 


0.30 — — 
1.14 0.84 0.00 
—0.57 —0.83 —0.01 
0.23 0.00 0.00 
0.26 —0.01 —0.01 
- 0.37 —O.12 0.07 
0.17 0.18 0.04 
0.48 —0.12 0.17 
0.08 0.07 0.11 
—0.04 —0.02 0.05 
0.55 0.10 0.03 





* Only the values of cloud type at 0000, 0300 and 0600 UTC are changed. 
Key — see Table IV. 





Table Vi. Typical values of thermal properties of road materials 





Source Substance Conductivity 


Heat capacity 
Wm' K"' . 


10°°Jm°K 





Farmer and 
Tonkinson 
(1989) 


Asphalt: wearing 
basecourse 
roadbase 
sub-base 
subsoil 
roadbed 
sub-base 
subsoil 


2.1 
2.1 
2.0 
1.9 
2.5 
2.4 
1.9 
2.5 





Concrete: 


1.40 
1.00 
1.00 
1.25 
1.30 
2.60 
1.25 
1.30 
1.30 
2.01 


Thornes 
(1984) 


Asphalt 
Concrete 


1.855 
1.675 


Oke 
(1987) 


Asphalt 
Concrete: 


0.75 
1.51 


1.940 


dense 2.11 


McCabe and 
Thompsen 
(1989) 


Weng et al. 
(1981) 


Individual B-road: 
laboratory 
field 

Multi-layer road: 
unfrozen 
dry ground 
wet ground 


Loam soil 70% saturated 


1.40 
0.97-1.48 


0.88-1.09 
0.63-0.95 
0.98-1.26 


1.58 


1.73-2.49 
1.59-2.40 


1.77-2.10 








Table Vil. 


Conductivity (k) (W m'' K"') and capacity (Cm) 
(10° J m’ K"') of road sub-layers for sensitivity tests 





Test 


Asphalt 
k Ce 


Concrete 


k 


Cm 


Soil (granite) 
k 


Cn 





Standard 


Test (a) 
Test (b) 
Test (c) 
Test (d) 
Test (e) 
Test (f) 


seuxses 
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Table Vill. 
Hill, 4 December 1987 


Results of sensitivity test for road thermal properties for Chapman’s 





Tests Output of tests 


SD 


Running difference 


DT(min) B Cc 





0.289 
0.269 
0.307 
0.278 
0.302 
0.289 
0.286 


Standard input 
Test (a) 
Test (b) 
Test (c) 
Test (d) 
Test (e) 
Test (f) 


0.00 
0.10 
0.00 
0.07 
—0.01 
0.04 
0.02 


—0.02 
0.02 
—0.01 
0.01 
0.00 
0.00 








Table IX. Results of sensitivity test for road thermal properties for Chapman’s 


Hill, 24 February 1988 





Tests Output of tests 


Bias SD 


Running difference 


DT(min) A B is 





—0.16 
—0.10 
—0.20 
—0.08 
—0.22 
—0.14 
—0.16 


0.573 
0.559 
0.590 
0.560 
0.595 
0.572 
0.567 


Standard input 
Test (a) 
Test (b) 
Test (c) 
Test (d) 
Test (e) 
Test (f) 


0.06 
—0.04 
0.08 
—0.06 
0.02 
0.00 


0.10 
—0.10 
0.11 
—0.16 
0.00 
—0.02 


0.30 
0.40 
0.20 
0.41 
0.14 
0.30 
0.28 





cloud amount and cloud type. Some error in the forecast 
and observation of cloud cover is usually unavoidable 
and is thus one of the major sources of model prediction 
error. 

The influence of dew-point temperature and wind 
speed are normally small. The impact of wind speed 
on model output is greater on ‘extreme’ days than on 
‘damped’ days. A 20% variation in the thermal 
conductivity and capacity of the layers of road 
construction and soil has a small impact on model 
predictions. 
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Appendix 

In spectral analysis, phase (P;.(/)) and coherence 
(C,(/)) are the most useful parameters for measuring 
the relationship between two time-series x;(t) and x;(1). 
Here / stands for frequency which is the reciprocal of the 
length of time (period) for x(t) to go through one 
complete cycle. The phase and coherence are calculated 
from the power spectral densities (f,;(/), fi.(/)), cospect- 
rum (c;,x(/)) and quadrature spectrum (qg;:(/)) by means 
of the relations 


cil) + gil) 


2 oe ata 
ll) = FD full) 


(Al) 


Pi Al) = — arctan(q;/)/ c/)). (A2) 
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The spectral densities, cospectrum and quadrature 
spectrum in equations (A1) and (A2) are estimated from 
cross-correlation functions (Koopmans 1974, Bath 
1974). 

The phase function is simply the difference of phase of 
the two time-series at frequency / or, more precisely, the 
phase lead of x;(#) over x(t). In practice it is expressed 
by a time difference in the form 


tl 
The coherence (coefficient of coherence) has the 
properties of the absolute value of a correlation at 
frequency /. In some sense it measures the degree of 
linear association between the time series x;(t) and x;(¢). 


The significance of coherence is assessed by the 
statistical variable (Koopmans 1974) 
v-)N)C2 
potas 
1~¢y 
where v is calculated from 


7 2n — m/2 


m 


Here n is asample number or size, m the maximum lag 
and C-. acritical value. F’ follows the F-distribution with 
2 and 2(v — 1) degrees of freedom. If 

Cu >C? 


the coherence is considered to be significant. 


5§51.510.42:551.511.61:628.53 





A scheme for estimating fluctuations in 
concentration of an emitted airborne pollutant 


F.B. Smith 
Meteorological Office, Bracknell 





Summary 





A simple practical scheme is presented for estimating the probability of exceeding a prescribed concentration, 
when averaged over a given time-interval, of a passive pollutant emitted from a finite-sized continuous source. A 
brief summary of the background theory and typical examples are given, with the option of obtaining the required 
answers simply using the diagrams presented or, somewhat more accurately, the analytical equations. The method 
has widespread applications in risk analyses wherever the danger of accidental release of hazardous materials to the 


atmosphere exists. 


1. Introduction 


Fig.l shows an exampie of a time-record of 
concentration in a plume at a point some distance 
downwind of a continuous point source. It shows that, 
as expected, the plume must be meandering under the 
action of eddies that are bigger than the width of the 
plume, so that the concentration is intermittently zero 
over short periods whereas at other times the concen- 
tration is other than zero, albeit quite variable. This 
complex behaviour can be very important for many 
hazardous pollutants. Consider three examples: 

(a) The pollutant may be inflammable within a band 

of concentrations, outside of which either there is 

insufficient material to maintain ignition or too little 
oxygen. Consideration of the mean concentration 
alone might suggest this dangerous band is restricted 
to a fairly short range downwind. Consideration of 
the distribution of concentrations, the fluctuations, 
would however indicate that inflammable concen- 


trations could occur over a much greater range of 
downwind distance. 

(b) The pollutant may be toxic (e.g. chlorine gas) and 
the effective toxicity may increase more rapidly than 
linearly with concentration c. For example, the 
toxicity of chlorine is believed to increase like c*”’. 
Thus, whilst consideration of the mean concen- 
tration field might suggest it is fairly safe to enter the 
plume at a given location, in reality, because of the 
fluctuations, it would be highly dangerous since one 
breath of chlorine at the concentration of one of 
the peaks could be fatal. The non-linearity further 
emphasizes the danger. 

(c) Some odours can only be perceived by the human 
nose above some threshold. The designers of an 
industrial plant that has some odorous pollutant as a 
waste product must take into account the likely 
fluctuations in concentration outside the boundaries 
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of the plant if they are to avoid a stream of complaints 
from the public. Rightly, they are not concerned 
whether the mean concentration can be smelt, but 
rather whether they can detect the peak concen- 
trations. 


Most progress has been made so far in specifying the 
statistics of the concentration fluctuations with time at a 
point. However, there is much interest in the spatial 
distribution of concentration fluctuations. For example, 
with an inflammable gas, if ignition occurs at some point 
in the plume where the concentration at the time is 
within the flammable range it is important to know what 
is the probability that the flame could travel back along 
zones of the plume within the range of flammability 
without any breaks back to the source, or to some 
location of particular concern (like a home). This is a 
very difficult problem, although Thomson (1990) has 
made progress in specifying the nature of spatial 
covariances of concentration. 


2. Available research 

It is only within the last 10 years or so that effective 
research has been directed at this problem, and only now 
that the beginnings of a practical scheme are emerging. 
This is the result of complex theory and experimental 
data. The scheme will no doubt be improved with time, 
but already it is able to give very useful guidance on the 
magnitude of the fluctuations and the probability of any 
concentration being realized at any one time at any 
downwind position. The background theory will not be 
given in any detail here. For further details see Sawford 
(1987), Sykes (1988), Thomson (1990) and Mylne and 
Mason (1991). 


3. The plume and the basic concentration 
fluctuation parameters 

Although the scheme is intended for use with plumes 
in the atmospheric boundary layer, it will assume 


Example of a record of concentration at a point downwind from a point source showing the nature of the concentration 


homogeneous isotropic turbulence with an exponential 
Lagrangian velocity correlation function with finite 
time-scale r,. These assumptions are the easiest to 
study, and the results provide a useful test case for 
verifying and developing ideas. Moreover, in the 
atmosphere, several stages of the plume’s development 
occur before the inhomogeneities become important. At 
large travel-times the inhomogeneities (or the non- 
stationariness) always become important so that the 
user must be wary of the results at long range. Initially 
we will also assume the source has sufficiently small 
dimensions that it can be treated as a point source. 

If o; is the standard deviation of single particles in the 
across-wind y direction then 


oi = 2A(exp(—T) — 1 +T) (1) 
where A = o-ri and T= t/r1. 
If or is the variance of the y component of the 


centroid of a pair of particles taken about the centreline, 
then 


of = 4A [exp T)—exp(-7)+2+ 


+3f1-e(-Z)}], 


If 204 is the variance of the y component of the 
separation of a pair of particles, then 


a = 44 [exp (-) - i> 


Ce ee Ae ] 
* 773 ( _ ( 2 } 
Equations 1-3 satisfy the identity 


2 _— 2 
04 + o: = 203. 
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The forms of these three variances are shown in Fig. 2. 

Thomson’s random walk studies have supported the 
following relationship between the concentration fluctu- 
ation and the mean concentration at any point on the 
downwind centreline: 


(0) = uaXo)'( 2 (5) 


2 
Oa 
h 1 + 2ex 722.) and where (0) 
= —_— — Tr =——... 
where yu € p 3VA Cc aweewt : 


From these relationships we can calculate uc(0)A/Q 
and a-/c(0) as functions of t/r., where o- is the root- 
mean-square concentration fluctuation on the centreline 
given by 


(6) 
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Figure 2. 
T=¢/T1. 


Forms of the three variances defined in the text with 


4. The concentration probability density 
function 

Observations made to date indicate that the probability 
of having a concentration c at any instant is given by the 
so-called clipped-normal probability density function 


(PDF) 


p = at, exp{-3(-7) } 


if c>0, or 
| +) 
sy 5{! ef (35) 


if c=0, indicating that there is usually a finite fraction 
of the time, called the intermittency, when the plume is 
off the centreline and c = 0 (note, not all papers define 
intermittency this way). A clipped-normal distribution 


is a Gaussian (or normal) distribution that is centred © 
off the origin, and the part to the right of the origin 
stands, but the part to the left is integrated and 
represented by a 6-function of the same magnitude 
placed at the origin to represent the probability of 
having zero concentration. This distribution satisfies the 
requirement that _, f “p(c)de = |. The two parameters o 
and y are determined from 


- freee B)eatg oo(-Z) 0 


2 —2 
Oct 


= $(I +¥)f1 + erf()}+ 


+e exp (- x) 


o 


(8) 


x 


erf(x) = = | exp(— s*)ds. 
0 


Values of erf(x) are given in Table I. 

The values of o/¢ and yy are determined as functions 
of o./¢, which in turn has already been determined as a 
function of 7. Having thus found y as a function of T, 
the magnitude of the intermittency / is obtained from 
the PDF: 


ie Bhs +) 
r=3fi ert (+ ; (9) 
These various functions are plotted in Fig. 3. The 


variation of o./¢ with T is like 7“ for T< 2 and like 7 
for T> 10. This can in fact be seen most easily in Fig. 2. 
a/é 
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Figure 3. Variation of the concentration fluctuation parameters 
with travel time on the averaged plume centreline assuming a clipped- 
normal probability distribution function. 
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Table I. 


The error function erf (x) extracted from Jahnke and Emde (1945) 








000 
125 
227 
286 
284 


113 
236 
335 
389 
380 


292 
117 
847 
480 
*019 


205 
039 
778 
421 
969 


194 


530 


427 
802 


468 
835 
103 130 
340 361 400 

23 39 54 69 83 


508 
868 


548 
900 
205 


61 73 84 95 *06 
63 72 80 88 96 
38 44 50 56 61 
91 95 89 *03 

28 31 34 37 39 


32 52 72 91 *09 
02 15 28 4l 53 
14 22 31 39 46 
86 91 97 *06 
31 35 38 41 La 


59 61 63 65 67 
76 78 79 80 81 
87 87 88 89 89 
25 29 33 37 41 
59 61 64 66 68 





676 
790 
869 
893 
847 


789 
900 
974 
992 
937 


013 
*118 
183 
*187 
*117 


716 
494 
*175 
761 
*254 


798 
566 
*238 
814 
*299 


959 
708 
*361 
918 
*385 


624 
961 


661 
991 
229 252 
438 456 

97 “a 


698 
*020 
275 
473 
*24 


768 
*076 
319 
507 
*49 


*048 
297 
490 
*37 


*16 *26 *36 *45 *55 
*04 “43 *18 *25 *32 
67 72 77 82 86 
*11 *15 *18 *22 a 
42 a5 47 49 51 


*26 *42 *58 "73 *88 
64 75 85 95 *05 
54 61 67 74 80 

*11 *15 *20 *24 *28 
47 50 52 55 57 


69 71 72 74 75 
82 83 84 85 86 
90 91 91 92 92 
44 48 51 54 56 
70 72 73 75 77 











If we suppose that there is some threshold value of c 
above which detrimental consequences occur then the 
probability P that c>c» (where cs is this threshold 
value) is given by 

i: 


od | (es_ 
P= 3 [ erf { Ya \e 
5. Example 

An example may help to see how the scheme can be 
operated. Suppose u=6 ms", r,=6s,0,=0.5ms", 
the source strength Q= 100gs' and the downwind 
distance x = 120m. What is the probability that the 
concentration c exceeds 200 mg m”™ there? 

From the given parameters, A=9, T= 3.333 and 
pA/Q=0.54. Therefore from equations (1), (2), (3), (5) 
and (6) or Figs 2 or 3: 


(10) 


(a) logio0.54¢ = —1.474. Therefore ¢ = 62.2057 

mgm. 

(b) oc/t= 1.2431. Therefore o. = 77.35 mgm”. 
Using equations (7), (8) or (9) or Fig. 3 we infer 

(c) y = 0.25106, 

(d) o/@= 1.8623. Therefore o = 115.84mgm°”, 

where, if Fig. 3 is used, the values of the parameters 

are found for T= 3.333. 


The intermittency then works out to be 0.4 and 
P(c > ce) = 7.28%. 


6. Concentrations off the centreline 


The profile of concentration is assumed to be normal 
in the crosswind direction with standard deviation o;: 


Q (- yt? 
2rrusi aes 20% 

Thomson’s random walk studies indicate that across 
the plume @ ” behaves like 


¢-= 


(11) 


(0) (E/(0)* 
where g = 203/02. gis approximately | out to 7= 2 and 
thereafter slowly increases asymptotically to 2. 0. is then 


- 
c= 


(12) 


obtained from o2 = c’—¢*. The other parameters can be 


found exactly as before on the centreline, with the caveat 
that the value of P({c,) should be the smaller of that 
found by the method just described, and that found by 
assuming that P(c) is proportional to the ‘mittency’ m 
(to coin a new word meaning the opposite of 
intermittency J, i.e. m= 1 — I, the fraction of the time 
spent within the plume). 


P(cs(y)) _ mW) 


P(ca(0)) — m(0) (13) 
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The value of m can either be found from the method 
described earlier for points on the centreline, or by 
m(y)/m(0) = c(y)/ (0). The two estimates are somewhat 
different and, although the first method is more 
accurate, the second is more consistent with the 
suggested method of finding P(cs). 

In our example consider a point 10m off_the 
centreline at x =i20 m. Using equation (11) we get c(10) 
= 19.2574mgm™ and using equations (1), (2), (3), 
(6) and (12) of = 42.6421, 04 = 18.9569, o: = 66.3274, 
q = 1.2858 and c*(10) = 2180.86, o{10) = 42.5442, and 
a{10)/c(10)= 2.2092. 


Now we infer -y from equations (7) and (8) since we 
have o-/¢. Alternatively Fig. 3 can be used for the + and 
the remaining parameters using an artificial value of 
T=1.49 corresponding to o-/¢ = 2.2092. Therefore 
y = —0.5739. It follows from equations (7) to (9) that 
o = 109.45, J=0.717, and P(c > 200) = 0.76%, almost 
one tenth of the value of P on the centreline. Note that 
the alternative way of estimating P would have given 
2.25%, so taking the minimum estimate, the first value 
stands, although the degree of uncertainty is rather high 
at these very low values of P. 


7. The effect of averaging time 

The fluctuations typically have a much smaller time- 
scale than rz. Consequently the peak concentrations 
may only persist for a relatively short time. The effect of 
a Dollutant, and in particular a toxic pollutant, may be 
governed by the average concentration in one breath 
which takes one or more seconds to draw into the lungs. 
It is expected, then, that the peak averaged concentrations 
may be significantly smaller than the instantaneous 
peak concentrations. Random walk simulations suggest 


that the time-scale 7; for fluctuations at a point is given 
by 


og 
tr = 0.43—. (14) 
u 
The effect of averaging over a time ¢, is then to reduce 
the concentration variance according to the equation 


2 (ts) ri te ta 
0) 7 a {exp (- + ie +h. 


(15) 

In the example we have been following, at 120 m, let 
t2= 1s. Then o, = 4.354 and rr= 0.29. It follows that 
ot.) = 48.9. From Fig. 3 we infer y = 1.0 and o= 61. 
Hence P=0.06%, which is clearly a very significant 
reduction as a result of time averaging. 

Thomson has looked at the spatial correlation 
function for the concentration values. If 


_ exe{x +A) 
, SS ee 


Gc 


then r=1 at A=0, and it falls away like 


A 2s 
1-4(4) 


up to about A = oa, whereafter it behaves like (o4/ A)” 
out to at least A = L. This slow tail, which is only evident 
in the longitudinal direction (parallel to the mean wind), 
he ascribes to inertial meandering of the plume. 


8. The effect of source size 
The effect of source size is rather complex from a 


theoretical standpoint. In practice however, its main 


consequence is to reduce the concentration variance out 
to T= 1. This is expected because the finite size must 
broaden the instantaneous plume and thereby slow 
down the penetration of ‘clean’ air into the body of the 
plume, the process which is responsible for generating 
internal fluctuations in concentration. Referring back to 
the section on the concentration fluctuation parameters 
off the centreline, the finite size (denoted by oo) modifies 
both the formulae for g and c’(0); 


(16) 


(oi + 00) 
(oz + 03)(04 + 00) 





(0) = peo)’ (17) 

The consequences are summarized in Fig. 4 where 
experimental data due to Fackrell and Robins (1982 ) 
are compared with Thomson’s random walk modelling 
results. As can be seen in Fig. 4, the larger oo is the 
smaller is the maximum variance and the further 


downwind it occurs. These points are further illustrated 
in Figs 5 and 6. 


9. Real situations in the boundary layer 

For a source near the ground the assumption of 
isotropic homogeneous turbulence does not appear very 
realistic. However, the high frequency turbulence 
responsible for generating the international spread of 
the instantaneous plume and its internal concentration 
fluctuations is virtually isotropic and homogeneous, 
especially in the earlier stages of the plume’s growth. 
Broadly speaking then, we can assume the plume 
develops as predicted by isotropic theory until o,/T is 
about |, after which it evolves in a self-similar manner; 
that is, it is as though t/7, becomes constant since, as 
the plume deepens, so the mean value of rz increases. 

In very stable conditions, a difficulty arises in 
specifying the sometimes large meanders of the plume 
due to intermittent gravity-driven motions as pools of 
cooler air slide down into hollows and valleys. Likewise 
some difficulties occur in very unstable conditions due 
to the complexity of the buoyancy driven motions. 
Increased research interest is now focused on concen- 
tration behaviour in light-wind situations, whether in 
very unstable or very stable conditions. 
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Figure 4. Variation of normalized concentration variance with travel time and with source size. A comparison of theory due to Thomson 
(1990) (solid symbols), and experimental data due to Fackrell and Robins (1982) (open symbols). 
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Figure 5. Variation of maximum concentration variance with source size, according to experimental data due to Fackrell and Robins (1982). 
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Figure 6. Variation with source size of the travel time at maximum concentration variance (Fackrell and Robins 1982). 


10. Fluctuations in concentration at long 
range 

So far we have considered the fluctuations in concen- 
tration within a plume from a continuous source out toa 
distance of perhaps | km. Fluctuation analysis also 
becomes very important at downwind distances much 
greater than this when the time over which the emission 


occurs is smaller than the travel time, or when the 
synoptic situation changes over a time that is smaller 
than the travel time. In both these situations it is 
probable that any receptor only experiences a small part 
of the possible range of fluctuations. The concentrations 
the receptor sees, and the risk these may imply, can only 
be described in statistical terms. 
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Classically, plumes of debris have always been 
considered in the same terms as plumes of smoke 
emanating from a factory chimney. It was recognized 
that such a plume contains a lot of structure (as we have 
already seen), blobs and wisps of high concentration 
surrounded by areas of lower or zero concentration. 
However, when the plume is averaged over many 
minutes, these irregularities are largely ironed out and a 
smoother distribution of concentration is achieved. The 
plume can then be described in terms of just a few simple 
parameters. Principal of these are (i) the position of the 
centroid, or centre of mass, of the plume material 
(considered as a function of distance downwind from 
the source), and (ii) the width of the plume (again a 
function of downwind distance), often expressed as the 
standard deviation of the concentration distribution 
about the centroid. In reality the second parameter may 
be two parameters — the horizontal width and the 
vertical depth. 

Such a description has been used for plumes out to 
much greater distances than the short distances over 
which chimney plumes are normally visible. For 
example the plume from the aluminium smelter at 
Mt. Isa in Australia has been tracked by instrumented 
aircraft out to 1800km in simple meteorological 
situations when the plume’s character can still be 
described in terms of these few basic parameters (Carras 
and Williams 1988). They have analysed these measure- 
ments in terms of the way the width varies with distance 
(see Fig. 7). 


500 





x July 1979 
Vv April 1983 
© August 1983 


1 Divi rcai 


Plume width (km) 








Litiiil 
5 10 50 100 
Plume age (hours) 





Figure 7. Log-log plot of the measured values of plume width (as a 
standard deviation) as a function of plume age at Mount Isa, 
Australia. 


However, experience in the atmosphere, as well as 
modelling studies, suggest that such simple plume ideas 
are invalid in most meteorological situations after some 
20 hours of travel because by that time the plume has 
grown so wide that synoptic-scale variations in the 
wind-field can gradually distort and deform the plume 
into shapes that cannot sensibly be described by the 
simple parameters described above. In the conditions 
outlined at the beginning of this section, it is as though 
we are forced to look at the complex structure of the 
plume in a similar way as we looked at the instantaneous 
structure of the plume at very short range in earlier 
sections of this paper. 

This approach has been supported by our experiences 
with the radioactive-debris plume from the Chernobyl 
Reactor accident in late April and early May 1986. The 
resulting depositions, and the consequences they had on 
the safety of food and on wildlife stimulated many 
organizations and governments to review the models 
they had available at that time. Most, if not all of these, 
were found wanting and this resulted in new, widespread 
modelling activities. A brief description will now be 
given of one of these new models, that developed at the 
Meteorological Office. 


11. The UK Nuclear Accident Model 

The Meteorological Office has developed a sophistic- 
ated random-walk or Monte-Carlo model in order to 
simulate, in both forecast and hindcast modes, the 
transport, decay and deposition of radioactive debris 
emitted from any site in Europe. It is given the acronym 
NAME (The Nuclear Accident Modelling Exercise). A 
global version of the model capable of treating 
emissions from any site in the world will be developed 
shortly. The plume is simulated by releasing a large 
number of ‘particles’ at the source (at hourly intervals, 
reflecting the emission profile in time and space) and 
allowing them to be transported by 3-dimensional winds 
taken from the output of the Meteorological Office’s 
weather forecasting models. When run in the hindcast 
mode, the ‘actual’ analysed winds are used from earlier 
forecast runs rather than forecast winds. A 10-day roll- 
over archive is used to store past winds for this purpose. 
A time-step of 15 minutes is used (an appropriate eddy 
turnover time in the convective boundary layer) and at 
each step a random perturbation is added to the 
horizontal displacements (and to the vertical displace- 
ments outside the boundary layer) to account for 
turbulent diffusion: 


Xi0n = xj + u(x) At + Ar 


where At is the time-step, r is arandom number from the 
standard normal (Gaussian) distribution and A is a 
coefficient chosen so that the implied turbulent energy 
remains quasi-constant in time. Within the boundary 
layer the particles are randomly reassigned in the 
vertical at the end of each time-step so that over a period 
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of time each particle will sample the mean wind at each 
level. Particles are released from different heights to 
allow for initial buoyant rise at the source. 

The model is multi-level, with a realistically evolving 
boundary layer, inferred at each hour from the 
forecasting model output. As its height changes, so 
particles are entrained and detrained from the layer into 
the layer above. 

A mass of pollutant or radioactiyity is associated with 
each particle at the outset, and concentrations in air are 
calculated by counting the particles in each grid-volume 
weighted by the mass they each then carry. If the initial 
emissions are only coarsely known, comparison of the 
calculated concentrations can be compared with any 
field measurements to refine the emissions and the mass 
each particle carries. 

The particle ‘mass’ is constantly being modified by 
radioactive decay, by chemical or physical transform- 
ations, by dry deposition and by wet deposition. The wet 
deposition in the model requires a good knowledge of 
the precipitation field, and this is obtained by inputting 
the results of a semi-independent programme which 
combines observations from surface rain-gauges, from 
weather radar (corrected for false inputs), from satellites 
and from the weather forecast model’s output on rain. 

In areas where weather radar output is available the 
resolution of the rainfall field is as good as 5 km. Even if 
the model’s implied concentrations of pollution in the 
air contain some errors, at least such a highly resolved 














Figure 8. Simulation of the Chernobyl plume at midday on 2 May 
1986 using the UK Nuclear Accident Model (NAME). In spite of the 
relatively few particles simulating the release the model does quite well 
in indicating the debris over southern Europe and over eastern 
England. 


rainfall field will provide valuable information on where 
heavy depositions of pollution are likely to have 
occurred; and this is vital to agricultural monitoring 
teams in situations like Chernobyl. An example of the 
output of the model is given in Fig. 8. It confirms the 
distorted and convoluted nature of the plume. 
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Figure 9. Variation of the mean concentration within plume 
material as a function of distance. The concentration is equated to the 
number of particles found in the grid-square (roughly 35 km square) 
when 200 particies are released per hour. 
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Figure 10. Variation of concentration mittency (fraction of time 
the concentration exceeded the given threshold) with distance for 
several specified threshold dosages. Units of dosage are related to the 
sum over time of the number of particles in a grid-square in each 
time-step. 
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12. Statistical description of a plume at 
long range 

The NAME model has been run on several randomly 
selected occasions out to 5 days with a continuous 
release of particles from a specified source located at a 
gridpoint near the centre of the model area which lies in 
Cumbria. The plumes grow over these 5 days and show 
all the distorting features evident in Fig.8. The 
following parameters have been determined for both 
concentrations in the atmosphere, for time-integrated 
concentrations (or dosages), and for depositions on the 
underlying surface: 

(a) the mittency of the plume for specified threshold 

values, i.e. what fraction of the time the concentration 

exceeded the given threshold, 

(b) the mean concentration or dosage or deposition 

at specified positions and at specified times after the 

start of the release, and 

(c) the PDF of the ‘observed’ concentrations or 

dosages about their respective means. 


This work is not yet finished, but it already shows a 
degree of order and meaningfulness in its results. For 
example the distributions in (c) above are very similar to 
those at short range discussed in earlier sections. It still 
has to be clarified as to whether the PDF is best 
represented by the clipped-normal PDF implied at short 
range, or whether a sample exponential distribution is 
better 


P(c) = 1- exp(—$). 


Other examples of these preliminary results are 
shown in Figs 9 and 10. 
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Notes and news 


Fourth Workshop on Operational 
Meteorology, Whistier, B.C. Canada, 
15-18 September 1992 
The workshop, sponsored by the Atmospheric Environ- 
ment Service of Environment Canada and the Canadian 
Meteorological and Oceanographic Society, will have a 
principal theme of ‘Forecasting in the Nineties’. 
The Program Committee wishes to solicit papers 
on the following topics: 
(a) The meteorological data explosion — the integra- 
tion and effective use of information in an operational 
setting. 
(b) Forecast techniques and conceptual models — 
their place in forecast decision making. 
(c) Climate services — how can they be used 
effectively. 
(d) Delivery — techniques and requirements to 
effectively deliver forecasts to the user. 


The workshop format will consist of laboratory 
sessions, submitted papers, invited papers, panel 
discussions, poster sessions and demonstrations. A brief 
introduction of each poster session will be made during 
an appropriate oral session. 

Titles and reviewers’ abstracts of 400-800 words 
should be sent to 


Neil McLennan 

Chairman Program Committee 
Atmospheric Environment Service 
Suite 200, 1200 West 73rd Avenue 
Vancouver B.C. 

Canada V6P 6H9 


Authors should indicate their preference for presenting 
their paper orally, in alaboratory or poster session, or as 
a demonstration. Preferences will be considered to the 
extent possible. Abstracts will be evaluated on their 
relevance to the theme as well as on quality. Papers not 
related to operational meteorology will not be accepted. 
The deadline for laboratory submissions is 1 October 
1991, and for all others is | February 1992. Authors will 
be notified regarding the acceptance of their abstracts 
and instructions on the format of their papers by 
1 November 1991 for laboratories and by | March 1992 
for other sessions. 

Complete camera-ready papers of not more than 
eight pages, including diagrams, must be received by the 
Program Chairman no later than 15 June 1992. A pre- 
print volume will be prepared and distributed to all 
registered workshop attendees. 

For additional information contact either Neil 
McLennan (Phone: 604-664-9073 Fax: 604-664-9066) or 
Gerard Neault (604-664-9052). 


Reviews 


The telemetry of hydrological data by 
satellite, by I.C. Strangeways. 209 mm X 295 mm, 
pp. v+60, illus. Wallingford, Institute of Hydrology, 
1990. Price £7.00. 


This short book, it is only 60 pages long, is a 
contribution to the Institute of Hydrology Report Series 
(No. 112). It is softback, and is intended to provide the 
basic information that a prospective user of satellite 
telemetry might need in setting up a communications 
system using facilities offered by both geostationary and 
polar-orbiting satellites. The author regards the inform- 
ation contained therein as a handbook intended as an 
introduction to satellite techniques and as an encourage- 
ment to use them. He is well placed to write such a 
document having been involved in developing automatic 
weather stations and using satellite telemetry for many 
years. 

The information is contained in six chapters covering 
the Meteosat system, the practical evaluation of a 
satellite telemetry system, commercial equipment and 
operating costs, the international situation and conclus- 
ions. There are eleven references, of which five are to 
work reported by the author, and this is an indication 
that in reading this book one should bear in mind that it 
represents a rather personal view of satellite telemetry. 
In spite of this, it is extremely useful in that it sets down 
much material which is spread widely through a range of 
literature. The text is clear and easy to ready, and most 
of the diagrams are clear and concise, although some of 
the data listings have not reproduced well (for example 
fig. 3.6). 

The emphasis throughout is on the Data Collection 
Platform (DCP) facility provided on the Meteosat series 
of satellites. The Argos polar-orbiting satellite communi- 
cations system is mentioned, but there is no mention of 
commercial VSAT systems, the Meteosat Data Dissemin- 
ation Mission (MDD) or the ESA CODE experiment 
based upon the Olympus satellite. Nevertheless the 
concentration on the systems most used by meteorologists 
and hydrologists serves to emphasize the practical 
nature of the book. 

Astrong argument in favour of using satellite communi- 
cations is presented, based upon the falling costs of the 
equipment needed and ease of maintenance, although it 
is assumed that the costs of any such services in the 
future will remain free. It is argued that there are a 
significant number of benefits in using satellite telemetry 
and these are listed, although it is recognized that the 
proportion of these benefits with hydrological applica- 
tions is not known. The reader is referred to ‘an 
authoritative guide’ as to the future of satellite telemetry 
in hydrology obtainable from ESA. It is a pity that this 
book was produced before the issue of WMO Report 
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TD No. 256 entitled ‘Use of and requirements for 
satellite imagery and data transmission for hydrological 
purposes’ by T. Kinosita (WMO Tech. Reports to the 
Commission for Hydrology No. 24, 17 pp, 1990). This 
publication takes the form of a report on a questionnaire. 
It notes that out of 34 countries who responded 27 did 
not use satellites for the transmission of hydrological 
data from data collection platforms, whereas 16 
countries said that they had plans to use satellites in this 
way as follows, 
Satellite type: geostationary lI 
polar orbiting 9 
communication 4. 


This information has to be reconciled with that con- 
tained in the review publication, which notes on page || 
that in the 1990s it is expected that 5000 DCPs out of a 
total of 12500 will be used for the collection of 
hydrological data. The implication is that most of the 
usage of DCPs for hydrology may be concentrated in a 
few countries. At the time the author wrote his book. 
only 400 stations were registered to use Meteosat and 
1800 were registered with the US Geological Survey. 
This total of 2200 was expected to increase to 12 550. 
This expected increase is supported by the WMO 
survey, although the sixfold increase expected by ESA 
may be too optimistic and a three- or fourfold increase 
might be more realistic. The sentiment expressed by the 
author that ‘anyone contemplating satellite telemetry 
for hydrology will probably use a geostationary satellite’ 
(page 13) is not borne out by the WMO survey. 

Whilst the majority of the information contained in 
this book is accurate, and there are useful descriptions of 
the Meteosat ground segment facilities and message 
dissemination procedures, there are some inaccuracies 
and omissions. For example, Eumetsat is not strictly 
funded by just the member countries of ESA but by the 
National Meteorological Services. There is no discussion 
of what causes system downtime, for example satellite 
decontamination exercises, eclipses, etc. Data distribution 
and charging policy issues are not discussed. Cost 
comparisons with UHF and microwave ground telemetry 
systems are not made recognizing that in the future the 
acquisition of sub-daily data may not be cost-effective 
using satellite systems if the service attracts a charge 
from the satellite operators. 

Unfortunately the material contained in this manual 
will become out of date fairly quickly. For example, 
some of the information on Meteosat Second Generation 
is already out of date and the section on commercial 
equipment likewise. In spite of the limitations noted in 
this review, the author has assembled a lot of very useful 
information on satellite telemetry. This manual is an 
important source of practical data, and should therefore 
be essential reading for all those contemplating 
acquiring hydrological and meteorological data via 
satellites. 

C.G. Collier 


The hurricane, by R.A. Pielke. 184 mm X 252 mm, 
pp. xi+228, illus. London, New York, Routledge, 1990. 
Price £50.00. ISBN 0 41503 705 0. 


The hurricane is one of meteorology’s most fascinating 
and awesome phenomena, and its destructive force is 
seldom appreciated by those living out of the tropics. It 
is also an excellent example of many important 
meteorological processes at work — convection and the 
release of latent heat are essential components of this 
heat engine, as are the exchanges of sensible heat and 
moisture at the surface. Hurricane motion, so important 
to forecast accurately if lives are to be saved, is often 
irregular and depends not only on the large-scale 
steering flow in which the system is embedded, but also 
on its interaction with troughs and jets in middle 
latitudes. Any book devoted to the subject has to 
consider extratropical as well as tropical meteorology 
and has, therefore, plenty of ground to cover. 

The author attempts to do just that in a non- 
mathematical framework which is aimed at the educated 
layman. The publisher’s hand-out describes it as a ‘well- 
illustrated and accessible book providing a state-of-the- 
art assessment of the climatology, forecasting and 
physical processes associated with tropical cyclones’. 
Well illustrated it certainly is. having some 190 pages of 
figures and diagrams. However, the choice of material is 
very strange, with over half the book taken up by maps 
of the tracks of every hurricane in the North Atlantic 
since 1871; of historical interest perhaps but barely the 
basis of a well-balanced survey of the subject. If it is to be 
heavily weighted towards the pictorial, which is not a 
bad idea for a subject like this, why not have more 
satellite pictures showing not just the textbook hurricane 
with an eye, but also the stages of development and the 
different structures that can occur? There are only three 
satellite pictures in the book. and all are of a large and 
fully mature system. The range of topics covered in the 
text is wide, but the explanation of the meteorology, 
being condensed into about 30 pages, is necessarily very 
sketchy. In fact the book gives the impression of being 
put together rather rapidly with insufficient thought 
given to clarity of presentation. To take just one 
example, the reader will be left wondering how the four 
pages devoted to diagrams showing the mean wind at 
850 mb ‘demonstrate clearly the validity’ of criteria for 
hurricane development, in the absence of any further 
explanation. This only goes to show that writing at a 
level for the educated iayman is really far harder than 
writing for the specialist; clarity is all important as is the 
ability to appreciate the reader’s viewpoint in order to 
explain what it is that he doesn’t understand. In this 
respect, | am afraid, the book fails. 

C.D. Hall 
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Satellite photographs — 17 May 1991 at 0600 and 1300 UTC 











Figure 1. Meteosat visible images of the Persian Gulf area on 17 May 1991 at (a) 0600 UTC and (b) 1300 UT 
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Meteorological satellite imagery of the Persian Gulf 
area continues to show the effects of the many oil fires 
that are still burning in Kuwait. The smoke plumes have 
been seen to extend for many hundreds of kilometres 
and one example of this can be seen in the images shown 
here. 

The two reprojected visible Meteosat images (Fig. 1) 
are products from the Meteorological Office's dedicated 
satellite-image processing system, Autosat-2. They 
show the type of imagery used within the Nowcasting 
Section of the 
Division to monitor routinely the extent of the smoke 
over the Gulf region. This task was initiated in January 
this year, when smoke was first apparent in satellite 


Short-range Forecasting Research 


imagery. 

The plume shown (formed by the smoke from the 500 
oil fires reported to be still burning in the area), 
extended southwards at a rate of 15 kn, consistent with 
the winds below 850 mb (5000 ft), as shown on the 
tephigram (Fig. 2). 

In the 1300 UTC image (Fig. I(b)), cumulus clouds 
can be seen within the smoke area. These have been 
termed ‘pyrocumulus’ and form typically during the day 
in this type of situation. One possible explanation for 
the formation of these clouds is that the upper layers of 


the smoke plume (having an albedo of 5-8%, as 
measured in recent field studies by the Meteorological 
Research Flight) absorb solar radiation, causing local 
heating of the surrounding air and hence strong ascent. 
In this situation, the unstable atmosphere allowed the 
developing clouds to penetrate high into the atmosphere. 
perhaps as high as 300 mb (30 000 ft), as illustrated on 
the tephigram and as suggested by their low brightness 
temperatures on the infra-red imagery (not shown here). 

These clouds are unlikely to provide a mechanism for 
the transport of significant quantities of particles into 
the stratosphere. Firstly, precipitation from these clouds 
will have washed out the particles (typically between 0. | 
and 10 wm) to the surface: secondly, the tropopause ts 
very much higher than 300 mb. 

The images show a light grey region covering a large 
area to the west of Kuwait. Its appearance and 
movement, combined with the 0600 UTC observations, 
(Fig. 3) indicate that it was a sandstorm, and was 
moving at a speed consistent with the observed near- 
surface winds. 

The observations at 0600 UTC within the region 
covered by the image, reported generally very low 
visibility, even down to zero at some stations, indicating 
the extent and severity of the airborne sand and smoke. 

G. Holpin and R. Bosworth 
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Figure 2. Tephigram for 17 May 1991 at 1200 UTC for King 
Khaled International Airport (location marked ‘K’ in Fig. 3), Riyadh, 
Saudi Arabia (thick lines) (not through plume). Also shown is a 
schematic showing typical plume depth, and a postulated ascent (thin 
lines) through the plume (using Kuwait surface data), constructed on 
the assumption that insolation causes the upper layers of the smoke to 
behave as if they were a surface, resulting in ascent (dotted line) and 
cloud formation. 





Figure 3. Middle East plotted surface observations on 17 May 1991 
at 0600 UTC, with plume area shown stippled. 
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GUIDE TO AUTHORS 


Content 


Articles on all aspects of meteorology are welcomed, particularly those which describe resu‘ts of research in 
applied meteorology or the development of practical forecasting techniques. 


Preparation and submission of articles 


Articles, which must be in English, should be typed, double-spaced with wide margins, on one side only of A4-size 
paper. Tables, references and figure captions should be typed separately. Spelling shouid conforin to the preferred 
spelling in the Concise Oxford Dictionary (latest edition). Articles prepared on floppy disk (Compucorp or 
IBM-compatible) can be labour-saving, but only a print-out should be submitted in the firsi instance. 

References should be made using the Harvard system (author) date) and full details should be given at the end of 
the text. If a document is unpublished, details must be given of the library where it may be seen. Documents which 
are not available to enquirers must not be referred to, except by ‘personal communication’. 

Tables should be numbered consecutively using roman numerals and provided with headings. 

Mathematical notation should be written with extreme care. Particular care stiould be iaken to differentiate 
between Greek letters and Roman letters for which they could be mistaken. Double subscripts and superscripts 
should be avoided, as they are difficult to typeset and read. Notation should be kept as simple as possible. Guidance 
is given in BS 1991: Part I: 1976, and Quantities, Units and Svmbols published by the Reyai Society. S! units. or 
units approved by the World Meteorological Organization, should be used. 

Articles for publication and all other communications for the Editor should be addressed io: The Chiet 
Executive, Meteorologicai Office, London Road, Bracknell, Berkshire RG12 2SZ and marked ‘For Meteorological 
Magazine’. 


Illustrations 


Diagrams must be drawn clearly, preferably in ink, and should not contain any unnecessary or irrelevant details. 
Explanatory text should not appear on the diagram itself but in the caption. Captions should be typed on a separate 
sheet of paper and should. as far as possible, explain the meanings of the diagrams without the reader having to refer 
to the text. The sequential numbering should correspond with the sequential referrals in ihe text. 

Sharp monochrome photographs on glossy paper are preferred: colour prints are acceptable but the use of colour 
is at the Editor's discretion. 


Copyright 
Authors should identify the holder of the copyright for their work when they first submit contributions. 
Free copies 


Three free copies of the magazine (one for a book review) are provided for authors of articles published ta "1 
Separate offprints for each article are not provided. 


Contributions: It is requested that all communications to the Editor and books for review be addressed to the Chief Executive. Meteorological 
Office. London Road. Bracknell. Berkshire RG12 2SZ, and marked *For Meteorological Magazine’. Contribuiors are asked to comply with the 
guidelines given in the Guide to authors which appears on the inside back cover. The responsibility for facts and cpinicns expressed in the signed 
articles and letters published in Meteorological Magazine rests with their respective authors. 


Subscriptions: Annual subscription £33.00 including postage: individual copies £3.00 including postage. Applications tor postal subscriptions 
should be made to HMSO, PO Box 276, London SW8 SDT; subscription enquiries 071—873 8499. 


Back numbers: Full-size reprints of Vols |—75 (1866—1940) are available from Johnson Reprint Co. Ltd, 24—28 Oval Road. London NWI 
7DX. Complete volumes of Meteorological Magazine commencing with volume 54 are available on microfilm from Uriversity Microtilms 
International. 18 Bedford Row. London WCIR 4EJ. Information on microfiche issues is available from Kraus Microfiche. Rte 100, Milwood 
NY 10546, USA 
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